Thermoelectric transport through strongly correlated quantum dots 
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The thermoelectric properties of strongly correlated quantum dots, described by a single level An- 
derson model coupled to conduction electron leads, is investigated using Wilson's numerical renor- 
malization group method. We calculate the electronic contribution, K c , to the thermal conductance, 
the thermopower, S, and the electrical conductance, G, of a quantum dot as a function of both tem- 
perature, T, and gate voltage, v g , for strong, intermediate and weak Coulomb correlations, U, on 
the dot. For strong correlations and in the Kondo regime, we find that the thermopower exhibits two 
sign changes, at temperatures Ti(v g ) and T2(v g ) with Ti < T2. We find that Ti > T p (v g ) » Tx(v g ), 
where T p (v g ) is the position of the Kondo induced peak in the thermopower, Tk(v 9 ) is the Kondo 
scale, and T2 = 0(T), where V is the level width. The loci of Ti(v s ) and T2(v g ) merge at a critical 
gate voltage v g = Vg(U/T), beyond which no sign change occurs at finite gate voltage (measured 
relative to mid- valley). We determine v g for different U/T finding that v g coincides, in each case, 
with entry into the mixed valence regime. No sign change is found outside the Kondo regime, or, 
for weak correlations U/T < 1, making such a sign change in S(T) a particularly sensitive signature 
of strong correlations and Kondo physics. The relevance of this to recent thermopower measure- 
ments of Kondo correlated quantum dots is discussed. The results for quantum dots are compared 
also to the relevant transport coefficients of dilute magnetic impurities in non-magnetic metals: the 
electronic contribution, n e , to the thermal conductivity, the thermopower, S, and the impurity con- 
tribution to the electrical resistivity, p. In the mixed valence and empty orbital regimes, we find, as 
a function of temperature, two peaks in K e as compared to a single peak in k c , and similarly, G(T) 
exhibits a finite temperature peak on entering the mixed valence regime, whereas such a pronounced 
peak is absent in p(T) even far into the empty orbital regime. We compare and contrast the figure 
of merit, power factor and the extent of violation of the Wiedemann- Franz law in quantum dots 
and dilute magnetic impurities. The extent of temperature scaling in the thermopower and thermal 
conductance of quantum dots in the Kondo regime is discussed. 

PACS numbers: 72.10.Fk, 72.15.Qm, 72.15.Jf, 73.63.Kv 
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I. INTRODUCTION 

Materials with potentially useful thermoelectric prop- 
erties are currently under intense theoretical and experi- 
mental investigation, mainly due to the prospect of appli- 
cations, e.g., for conversion of waste heat into electricity 
in thermoelectric generators, for applications to refrig- 
eration, or for on-chip cooling and energy efficiency in 
microelectronics applications^— . Apart from possible 
applications, thermoelectric materials can also serve as 
an interesting testing ground for theoretical approaches 
to electrical and thermal transport in solids^—. As the 
scale of the individual components in semiconducting de- 
vices is approaching the nano-size, a description of ther- 
mal transport through quantum dots is also attracting a 
lot of experimental and theoretical attention 2 ^— . 

In this paper we address the thermoelectric properties 
of a nanoscale size quantum dot exhibiting the Kondo 
effect, which we describe in terms of a single level Ander- 
son impurity model with two conduction electron leads 
at fixed chemical potentials. The quantum dots that we 
consider have sizes of 10 — 100 nm and can be tuned 
from the Kondo to the mixed valence and empty or- 
bital regimes by a gate voltage^—. Short segments of 
carbon nanotubes2£ connected to leads exhibit similar 
physics, so our results could also be of relevance to such 



systems. Very recent experiments on nanoscale quan- 
tum dot s 20 ' 28 are beginning to probe the effect of Kondo 
correlations on the thermopower, although as we shall 
argue in the conclusions a quantitative comparison with 
theory is still some way off. The thermoelectric proper- 
ties of dilute magnetic impurities in non-magnetic met- 
als, such as Ce^L/ai-a; AI3 and Ce^Lai-^Bg, are closely 
related to those of quantum dots (see Sec. [TTJ) so we dis- 
cuss these here also. Understanding the thermoelectric 
properties of magnetic impurities is also a useful starting 
point for understanding those of heavy fermions within 
the dynamical mean field theory approach^ although, in 
these systems, crystal field effects and non-resonant chan- 
nels play a crucial role for the thermopower, and need to 
be taken into account for a quantitative comparison to 
experiment^—. 

The approach that we use in this paper, Wilson's nu- 
merical renormalization group (NRG) method^ - — , gives 
reliable results for transport properties in all parameter 
and temperature regimes of interest^. The present cal- 
culations were carried out for the Anderson model with 
finite Coulomb repulsion, as is appropriate for nanoscale 
size quantum dots. We implemented recent develop- 
ments in the calculation of dynamical quantities within 
the NRG, including the use of the self-energy^ and 
the full density matrix (FDM) generalization^ (see also 
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Rcf. 34.35) of the reduced density matrix approach^ 
within the complete basis set of eliminated states^. In 
particular, the FDM approach allows calculations of dy- 
namical properties at all excitation energies w relative 
to the temperature T, thereby simplifying the calcula- 
tion of transport properties which require knowledge of 
excitations, oj, above and below the temperature^. 

The outline of the paper is as follows. In Sec. ITlllnll wc 
describe the Anderson impurity model for quantum dots 
and dilute magnetic impurities and we specify the rele- 
vant transport quantities that we calculate for these two 
different physical realizations of the model. The NRG 
method used in this paper is described in Sec. QV] to- 
gether with results for occupancies which we use to de- 
fine Kondo, mixed valence and empty orbital regimes in 
the strong correlation limit. Sec. |V] presents the temper- 
ature dependent transport properties of quantum dots 
and Sec. I VII compares these to the corresponding quan- 
tities for dilute magnetic impurities. Results for the fig- 
ure of merit, power factor and Lorenz number ratios for 
quantum dot and magnetic impurity systems are pre- 
sented in Sec. I VIII Sec. IVIIII investigates the extent 
to which universal scaling functions apply to the ther- 
mopower and thermal conductance of quantum dots in 
the Kondo regime. In Sec. [IX] we present our results 
for the gate voltage (local level) dependence of transport 
quantities for quantum dots (magnetic impurities). Con- 
clusions and a discussion of the relevance of our results 
to recent experiments on nanoscale size quantum dots is 
presented in Sec. |Xj Appendix [A] discusses the reduc- 
tion of the two-lead Anderson model to a single chan- 
nel model, Appendix IBl contains some additional results 
for moderately and weakly correlated quantum dots, and 
Appendix [C] provides details of the FDM approach^ and 
an alternative detailed derivation of the FDM expression 
for local Green's functions, which we have used to ob- 
tain the results in this paper. Finally Appendix [Pi gives 
an outline of the derivation of thermopower and thermal 
conductance for quantum dots. 



states \a). We assume a flat density of states of magni- 
tude Np = 1/2D per spin channel for both leads, where 
D = 1 is the half bandwidth of each lead. The single- 
particle broadening (half-width at half-maximum) of the 
d-level is given by F = Fl + Fr, where r^fl = i^Npt 2 L R 
are the contributions to the broadening from the left and 
right leads. In this paper, we follow the convention used 
in quantum dot work and use as unit of energy not T, 
but the full-width at half-maximum F = 2F. 

Since the d-state of the quantum dot in ([1]) only couples 
to the even combination a eka ~ ^h^Lko + ^rCrucj of the 
lead electron states, one can show (see Appendix [A)) that, 
to a very good approximation, the above model can be 
reduced to the following single-channel Anderson model 

H = e ekcr a\ k(j a eka + ^ e d d\d a + U n d fn d i 

ka <y 

+ tY J {a\ ha d, + h.c.), (2) 

fc<7 

where the tunneling amplitude t is given by t 2 = t\ + t R , 
so that the hybridization strength of the dot to the leads 
is given by T = Fl + Fr. This is also the appropriate 
model for describing dilute magnetic impurities in non- 
magnetic metals^. In fact, for both systems (see below 
and Appendix [AJ the calculation of the linear transport 
properties reduces to the calculation of the equilibrium 
d-levcl spectral density of the single-channel model 

A(u) = --Ira[G d(T (oj + i6)], (3) 

where G d(T {uj+i8) = {{d a ;d\)) is the Fourier transform of 
the retarded d-level Green function of ©. Hence, all re- 
sults in this paper, including those for dilute magnetic 
impurities, are obtained by solving the single-channel 
model © using the NRG (as explained in Sec. IIVI) to 
obtain A(u,T). 



II. MODEL 

A nanoscale quantum dot is described by the single 
level Anderson impurity model with two conduction elec- 
tron leads 

H = e a kacl ka c akcr + y~] e d d^d^ + U n d ^n di 

aka 

Here, e ak( y is the kinetic energy of conduction electrons 
with wavenumber k and spin a in lead a = {L, R), e d is 
the local level energy, U is the Coulomb repulsion on the 
dot and t a is the tunnel matrix element of the dot level 
to conduction electron states in lead a = (L,R). The 
operators c^ fco .(c a fe CT ) create (destroy) conduction elec- 
tron states \aka) and d\ r {d <7 ) create (destroy) local d-level 



III. TRANSPORT QUANTITIES 

A. Quantum dots 

Thermoelectric transport through the quantum dot ([1]) 
is calculated for a steady state situation in which a small 
external bias voltage, SV = Vl — Vr, and a small temper- 
ature gradient ST is applied between the left and right 
leads. Left and right leads are then at different chemical 
potentials /xl and /x^, and temperatures Tp and Tr, with 
e5V = n L -ii R and ST = T L - T R . We follow the ap- 
proach for deriving the electrical conductance, G(T), the 
thermal conductance, K e (T), and thermoelectric power, 
S(T), through an interacting quantum dot^r— using the 
non-equilibrium Green's function formalism. For com- 
pletenes, an outline of this derivatio n 21 ! 22 can be found 
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in Appendix [D] The final expressions are given by 
G{T) = e 2 I Q {T) 
S(T) = 



1 h(T) 



K e {T) 



1 

T 



\e\TI (T) 
h{T) 



Io(T) 



where /„, n = 0,1,2 are the transport integrals 
UT) = \ J <W l 7»(-|JJ). 



(4) 
(0) 



(7) 



Here, e denotes the magnitude of the electronic charge 
and h denotes Planck's constant. The quantity is 
related to the spectral density A{uj) via 



T(w) = 4tt . L ? A{u). 
T L +T R 

At T = 0, the conductance acquires the value 
2e 2 



G(0) 



;t -r(o) 

2e 2 4f L f fl 
& (f L +f fl ) : 



■sin 2 (7rn d /2) 



(8) 

(9) 
(10) 



where n<2 is the occupancy of the dot and we have used 
the Friedel sum rule, 



A{lo = 0,T = 0) 



1 



n(T, 



■T R ) 



sin 2 (7rn rf /2). 



For integer occupation, n& = 1, and equal coupling to 
the leads, l\ = Tr, the conductance reaches the unitary 
value 2e 2 /h, which we henceforth denote by Gq. 



B. Dilute magnetic impurities 

It is of interest to compare the transport properties of a 
quantum dot, with the corresponding quantities for elec- 
trons scattering from a dilute concentration, n% <C 1, of 
magnetic impurities in a clean host metal with constant 
density of states Np per spin. As for quantum dots, 
the relevant model for such dilute magnetic impurities is 
the single channel Anderson model ^ with hybridization 
strength T. In order to obtain the thermopower, S(T), 
the thermal conductivity, K e (T), and resistivity, p(T) (or 
conductivity a = 1/p) for such a dilute concentration of 
magnetic impurities we use the Kubo formalism, see Ap- 
pendix A of Ref. UH for the details, and find for these 
quantities 



P(T) 
S(T) 

Ke(T) 



1 



l Mo{T) 



1 

T 



1 Mi(T) 
|e|rM (T) 

M 2 (T) 



M 2 (T) 



M (T) 



(11) 
(12) 

(13) 



The transport integrals M n ,n = 0,1,2, appearing here, 
are now defined by 



M n (T) 



du> u n T{w,T){-^-), (14) 



where t(u), T) is the transport time of electrons, which is 
given in terms of the impurity spectral density A{io,T) 
by 



1 



t{w,T) 



^2TA(uj,T). 

iv p 



(15) 



In order to compare impurity transport properties with 
those of quantum dots, using the same units, we shall 
use rescaled quantities, e.g. for quantum dots G/Gq and 
K e /Go and for impurities p/ po and n e po where 



po = Iniji^Npe 1 



(16) 



is the unitary resistivity of electrons scattering from a 
dilute concentration of magnetic impurities. 

While the physics governing the transport properties 
of electrons scattering from dilute magnetic impurities, 
described by ([2]), is expected to be similar to that govern- 
ing the transport properties of electrons through quan- 
tum dots (also described by (|5J)), differences are also ex- 
pected, particularly for the respective thermopowers or 
the thermal conductance (conductivity), for the follow- 
ing reason: the transport expressions for quantum dots 
arise from integrals /„, n = 0, 1, 2 which involve the n'th 
moments of A(oj,T) convoluted with the derivative of a 
Fermi function, whereas those for magnetic impurities 
arise from n'th moments of 1/A{uj,T) convoluted with 
the same derivative. At low temperatures, a Sommer- 
feld expansion for 1\ and M\ results in different signs 
for the thermopower in the two different situations, since 
derivatives of A and 1/A have opposite signs. On the 
other hand, at higher temperatures, moments of A and 
1/A are determining factors for transport. Particularly 
the moments I\ (Mi), entering the thermopower, and 
h (M 2 ), entering the thermal conductance (conductiv- 
ity), probe differences in the behavior or A(u,T) and 
1/A(uj,T) at high temperature. Consequently we expect 
significant quantitative differences for the thermopower 
and thermal conductance (conductivity) of quantum dots 
and dilute magnetic impurities at high temperatures. We 
discuss these differences in Sec. [VI] 



IV. NRG APPROACH 



A. NRG and dynamical quantities 



We calculate the spectral function A(uj,T) and the 
transport properties of quantum dots, by using the NRG 
approach^—. This method is numerically exact and can 
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FIG. 1: Temperature dependence of the occupancy for dif- 
ferent values of the gate voltages v g > in the Kondo 
(solid lines), mixed valence (dashed lines) and empty or- 
bital (dotted lines) regimes. We define these regimes by 
\n d (T = 0) - 1| < 0.25, \n d (T = 0) - 0.50] < 0.25 and 
\n d (T = 0)| < 0.25, respectively. 



2 
1.5 
1 



0.5 

0.8 
>0.4 



— ^— -L' 1 1 


(a) - 

U/T=8 








1,1, 


— t/t=o.ooi 
t/t=o.i 


i - — 




t/t=i 

— t/t=io 

— t/t=ioo 


(b) 


s \ 




u/r=8 : 


1,1, 


1 ; h ^-4^_- 




8 -6 -4 


-2 2 


4 6 I 



FIG. 2: Dependence on gate voltage, v 9 = (e^ + U/2)/T, 
of, (a), the occupancy of the quantum dot, rid , and, (b), the 
double occupancy, D — {ndfridj,}, at several temperatures. 



be used to calculate both static thermodynamic proper- 
ties as well as finite temperature dynamic and transport 
properties^. In brief, the NRG procedur o 29 ' 30 consists of 
the following steps, (i), a logarithmic mesh of e£ = DA~ n 
is introduced about the Fermi level ef = 0, and, (ii), a 
unitary transformation of the a e ka in ([2|) is performed 
such that /oo- = ^ k a e ka is the first operator in a new 
basis, f ntT , n = 0, 1, ... , which tridiagonalizes H c = 



E 



He EuEn=0 ^"(/n+lir/ni' + 

h.c.)., where the hoppings t n rj A~™/ 2 for a flat con- 
duction band 3 ^. The Hamiltonian with the above 
discretized form of the kinetic energy is now iterativcly 
diagonalized by defining a sequence of finite size Hamil- 



tonians 

H m = J~] eg d\d a + U ridfridi + 1 sXfogdg + h.c.) 

tr cr 
m— 1 



Yl tnifl+lafn* + h.C.) 



(17) 



n— 0,a 



for m > up to a maximum chain length N. For each to, 
this yields the excitations E™ and many body eigenstates 
\pm) of H m at a corresponding set of energy scales ui m 
defined by the smallest scale in H m . uj m = t m ~ A ~ . 
Since the number of states grows as 4 m , for m > 6 only 
the lowest 600 or so states of H m are retained. These are 
used as a basis for constructing H m+ i. For to > 6, both 
the retained and eliminated (high energy) states of H m , 
together with the corresponding eigenvalues, are stored. 
This information is subsequently used to evaluate the 
spectral function A(uj,T) within the FDM approach^ 3 - 
described in Appendix [Cj This evaluation makes use of, 
(i), the completeness of eliminated states 3 ^, allowing a 
multiple-shell evaluation of Green's functions^, avoid- 
ing double counting of excitations, and, (ii), the reduced 
density matrix approach to Greens functions, introduced 
to the NRG by Hofstetter— . In addition, we calculate 
the spectral function via the correlation part of the self- 
energy S(w,T) following Bulla, Pruschke and Hewson^ 2 -, 
via 



A(u,T) 
E(w,T) 



= U 



1 

- — Im 

7T 



1 



e d + iT - 



E(w,T) 
F a (w,T) 

Gdcr(iO,T) 



(18) 



Since the FDM entering the definition of the Green's 
functions, sec Appendix [Cl contains the complete spec- 
trum from all NRG iterations, asymptotically high and 
low temperatures can be investigated more easily than 
within previous approaches^, which involved at a given 
temperature T, choosing an appropriate energy shell to 
extract A(u),T). In addition, the regime ui <T, which 
was problematical in previous approaches, can now be ad- 
dressed, since contributions from all excitations (for all 
energy shells) are taken into account in the expression 
for the Green's function within the FDM approach. 



B. Calculations 

The calculations reported here have been carried out 
for a discretization parameter A = 1.75, retaining 660 
states per NRG iteration and a hybridization strength 
f = 0.01 (in units of the half-bandwidth D = 1). The 
maximum chain length diagonalized was N = 68. We 
use the full width r = 2r = 0.02 as our energy unit 
throughout. Results for a wide range of temperatures 
from T/r < 1 to T/T > 1 were obtained to fully char- 
acterize the transport properties of quantum dots and 
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dilute magnetic impurities. We note that, in practice, 
the regime T/Y 3> 1 is probably not accessible in experi- 
ment due to other effects which become important at high 
temperature, and which we do not take into account, e.g. 
phonons, multiple levels, crystal field states etc. Calcu- 
lations for strong (U/Y = 8,6), moderate (U/Y = 3) and 
weak (U /Y = 1) correlations were carried out for a range 
of dimcnsionless gate voltages, v g , defined by 

v 3 = £d = °- 25 «> n=±l,±2,..., ±32. 

With this definition, the gate voltage for mid-valley oc- 
curs at v g = for all U. Due to particle-hole symmetry, 
calculations were carried out for v g > 0, with those for 
v g < being obtained via a particle-hole transformation. 
This results in v g — > — v s , occupancy rid — > 1 — n d , dou- 
ble occupancy D = (rid^ridx) — > 1 — rid + D, thermopowcr 
S — > —S, with G and K c remaining unchanged. The 
behavior of S, G and K c under v g — > — v s , follows from 
their definition and the behavior of the spectral function, 
A(u,T) ->■ A(-u),T)) under v g -> -v 9 . 

Fig. [T] shows the temperature dependence of the dot 
level occupancy, rid(T), for gate voltages in the Kondo, 
mixed valence and empty orbital regimes, for U/Y = 8. 
For U/Y ^> 1 we use the occupancy at T = to de- 
lineate between the different regimes. Specifically, the 
Kondo regime is defined by gate voltages around mid- 
valley (v s = 0) with \rid{T = 0) — 1| < 0.25 (see caption 
to Fig. [T]). Similarly, the mixed valence and empty or- 
bital regimes are defined by gate voltages corresponding 
to \n d (T = 0) - 0.5| < 0.25 and n d (T = 0) < 0.25, 
respectively (see Fig.[T]). In the Kondo regime, a charac- 
teristic low temperature scale, the Kondo scale Tk, can 
be defined via^ 2 - 

Tk - U^e^ 1 = r/f e^-^)/*, (19) 

where u = U/Y. The mid- valley Kondo scales for u = 3, 6 
and 8, are T K /Y = 8.2 x 10~ 2 , 1.0 x 10~ 2 and 2.64 x 10~ 3 , 
respectively. 

Within the FDM approach, the thermodynamic value 
of the dot occupancy rid = Tr^n^], where p is the FDM 
defined in Appendix [C] and the value obtained from 
the spectral sum rule n d = ^ a J^°° — -^lm[F a (uj, T)]duj, 
with F a defined in (|18[) . are identical by construction, as 
we also verified numerically. Fig.[2]shows the gate voltage 
dependence of the dot occupancy (and for completeness 
also the double occupancy D = (n d ^ndi)) at a number of 
temperatures for the strong correlation case U/Y = 8. 

C. Physics of Kondo, mixed valence and empty 
orbital regimes 

Before presenting the results, a few words are in order 
concerning the physical significance of the Kondo, mixed 
valence and empty orbital regimes for strong Coulomb 



correlations (U/Y 3> 1) on the dot (for more detailed in- 
formation we refer the reader to Ref. [42h . The Kondo 
regime, nd ~ 1, corresponds to the formation of a lo- 
calized spin on the dot at intermediate temperatures 
(Tk < T <gc T). In this temperature range, physical 
properties exhibit logarithmic temperature dependences, 
the hall-mark of the Kondo effect. At T < T K the lo- 
calized spin is quenched by the lead electrons, resulting 
in a many-body singlet at T = and a narrow Kondo 
resonance (of width Tk) in the dot spectral density at 
the Fermi level. Physical properties are characterized by 
spin fluctuations on scales Tk < T <C Y, charge fluc- 
tuations on a scale T « Y and renormalized Fermi liq- 
uid excitations at T <C Tk- The dot spectral density is 
well understood^: it has a three peaked structure, with 
single-particle charge excitations at e d and Sd + U and 
a temperature dependent Kondo resonance at the Fermi 
level. The mixed valence regime, corresponds to gate 
voltages such that the level £d is within Y of the Fermi 
level. The charge on the dot fluctuates between n d = 
and rid = 1 resulting in an average charge n d ~ 0.5. The 
physics is governed by quantum mechanical charge fluc- 
tuations on a scale set by Y. The empty orbital regime, 
corresponds to rid ~ and Sd/Y 3> 1. Physical proper- 
ties are dominated by charge fluctuations, primarily via 
thermal activation (with an activation energy e d ). Even 
though U/Y 3> 1, the physics in this regime corresponds 
to that of a non-interacting resonant level model with a 
resonant level of width Y at energy Ed > 0. 



V. TEMPERATURE DEPENDENCE OF 
TRANSPORT PROPERTIES OF QUANTUM 
DOTS 

The temperature dependence of transport properties 
of a quantum dot described by the model ((T|) is shown in 
Fig. [3] for several values of the gate voltage, ranging from 
the Kondo regime (Fig. [3k,d,g), to the mixed valence 
(Fig. [3b,e,h) and empty orbital (Fig. [3b,f,i) regimes and 
for strong Coulomb correlations on the dot (U/Y = 8). 
Moderate to weak correlations are described briefly in 
Sec. IV Dl and Appendix 151 Depending on the regime, the 
transport properties exhibit different characteristic tem- 
perature dependences, which wc describe in detail below 
for each transport property in turn. Here, and in sev- 
eral other figures in the paper, we use arrows to indicate 
the evolution, with increasing gate voltage v 9 > about 
mid- valley (v g = 0), of the various transport properties. 

A. Electrical conductance: G(T) 

The general trends in the electrical conductance 
G(T) of Kondo correlated quantum dots are well 
understoo d 15 ' 43 " —: in short, as T — > 0, the conductance 
approaches a maximum value (see Fig. [3^), indicating 
that the quantum dot appears "transparent" to electrons 
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Kondo Mixed valence Empty orbital 




T/r T/r t/t 

FIG. 3: (a-c) The normalized electrical conductance, G/Go, where Go = 2e 2 /h, (d-f), the thermopower, S, in units of fes/e, 
and, (g-i), the normalized electronic contribution to the thermal conductance, K c /Go, multiplied by a factor 10 2 for clarity of 
presentation, as function of T/T, in the strongly correlated regime U/T = 8 for a range of gate voltages, v 9 = (ea + U/2)/T, 
in the Kondo (first column), mixed valence (second column) and empty orbital (third column) regimes. The range of v 9 is 
indicated in the top panels for each regime and the increment used was 0.25. Arrows indicate the evolution of the transport 
quantities with increasing v s . The inset to (d) compares the Kondo scale Tk in the Kondo regime with the peak position of 
the low temperature peak, T v , in S below the first sign change at Ti as a function of v g > 0. In the Kondo regime Tk ~ T p 
and for gate voltages approaching the mixed valence regime, the two scales deviate, as expected. The inset for K e in (g-h) 
shows the crossing point at T/T ~ 0.6 in more detail (vertical dotted line) and the evolution of the two-peaked structure for 
gate voltages approaching the mixed-valence regime. 



tunneling through it, and a logarithmic behavior around 
Tk marks the crossover from the weakly coupled regime 
at T > Tk to the strongly coupled regime at T C Tk- 
An issue, less discussed in the literature, which we point 
out here, is the appearance of a finite temperature peak 
in the conductance, G(T)/G$, on entering the mixed va- 
lence regime (see Fig. \3]p). This feature becomes par- 
ticularly pronounced in the empty orbital regime (see 
Fig. [3J;). This effect has been observed in experiments 
on lateral quantum dot s 23 ' 26 and a comparison to theo- 
retical calculations shows good agreemen t 48 ' 49 (see also 
the discussion of the resistivity of dilute magnetic impu- 
rities in Sec. IVI[) . 

B. Thermopower: S(T) 

The thermopower exhibits a particularly interesting 
temperature dependence in the Kondo regime, Fig. [5J1, 
with two sign changes at T = Ti(v s ) and T = T 2 (v s ), 
and, correspondingly, three extrema at T = T p , T w 
0.6 - 0.8r and T fts 6T. The detailed behavior of T x 
and T2 as a function of gate voltage will be described 
below; here, it suffices to note that neither T\ nor T2 



are low energy scales, and T2 is typically on a scale of 
order T (see Fig. [SJi and Fig. 2k below). The low tem- 
perature "Kondo" peak in S(T) at T = T p is found to 
scale with Tk (as defined in Eq. (jX9j)) , as shown in the 
inset to Fig. [5Ji. Thus, in contrast to T\ and T2, T p can 
be considered a low energy scale in the Kondo regime. 
The central positive peak in S(T) first grows with posi- 
tive magnitude on moving away from the Kondo regime 
(Fig. and then decreases in magnitude in Fig. \3jp on 
entering the mixed valence regime. Simultaneously, the 
"Kondo" peak in S(T) acquires a large negative value 
while merging with the high energy (negative) peak at 
T ~ 6r on entering the mixed valence regime (Fig. [3^). 
Well into the mixed valence regime, the thermopower 
exhibits a single negative peak on a scale T with a dis- 
tinct shoulder at higher temperatures due to the peak 
at T ~ 6r. This picture continues to hold in the empty 
orbital regime (see Fig.|3f), with the shoulder at T w 6r 
having almost disappeared. The thermopower remains 
negative for all gate voltages v g > in this regime. 

The above behavior in the temperature dependence of 
the thermopower in the Kondo regime is explained in 
terms of the structure of the single-particle excitations in 
A(u,T). At low temperatures, a Sommcrfcld expansion 
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for S(T) gives! 5 . 



Quantum dot 



Impurity 



1 



OA, 



S{T) -- k f\\ kBT A^ 



1=0: 



(20) 



showing that the sign of the thermopower depends on 
the slope of the spectral density at the Fermi level. For 
T <C Tjc and v g > 0, the Kondo resonance lies above 
the Fermi level, so the slope of the spectral density at 
the Fermi level is positive, resulting in a negative ther- 
mopower. This remains true on further increasing the 
temperature, but as shown in Ref. [TBI , eventually the 
Kondo resonance is suppressed at T > Tk resulting in 
a negative slope of the spectral density at w = for 
v g > (with the opposite being true for v g < 0). Con- 
sequently, the thermopower changes sign at the temper- 
ature T\ which roughly corresponds to the temperature 
at which the Kondo resonance vanishes. At T ^s> Tk, the 
determining factor for the sign of the thermopower is no 
longer the slope of the spectral function at oj = 0, but the 
number of states available below or above the Fermi level. 
These determine the overall sign of the transport integral 
Ii in the expression for the thermopower in Eq. ([5]). For 
v g > 0, there are rid/2 < 0.5 states below the Fermi level 
and 1 — rid/2 > 0.5 states above the Fermi level. Con- 
sequently, the integral of —(df/doj)uiA(u;,T) for uj > 
is greater than its counterpart for 10 < 0, so I\ > and 
the thermopower is again negative at T > Tk- This 
occurs at T = T2, which is found to be of order T (see 
below) . Due to the factor 1 jT coming from the derivative 
of the Fermi function in I\ , the negative thermopower at 
T > T-2 acquires a maximum negative value and then 
decreases as 1/T at T 3> T, exhibiting no further sign 
changes, as confirmed also numerically. We note that, 
away from half- filling (v fl = 0), the modified second order 
perturbation in U approac h 50 ' 51 gives an incorrect sign 
for the slope of the spectral density at the Fermi level 
in the Kondo regime. This results in a wrong sign for 
the thermopower at T < Tk in the Kondo regimes* com- 
pared to our NRG calculations (which agree with those 
of Ref. I52T ) . Approximate approaches using an infinite U 
Anderson mode l 53 ' 54 could also not access the low tem- 
perature Kondo regime. 

The sign changes of the thermopower of strongly corre- 
lated quantum dots at the temperatures T± and T2 in the 
Kondo regime, are particularly interesting. They provide 
a "smoking gun" signature for Kondo behavior in quan- 
tum dots, and could be used in future experiments as sen- 
sitive probes of strong correlations and Kondo physics. It 
is therefore interesting to give a detailed characterization 
of the dependence of T\ and T2 on gate voltage and in- 
teraction strength U/T. We show in Fig. [4^,-b the loci 
of T\ and T2 as a function of v g > for quantum dots 
and magnetic impurities for three interaction strengths. 
Although S(T) vanishes at v g = 0, T\ and T2 have finite 
limiting values there. These are difficult to determine 
numerically due to the vanishingly small thermopower in 
this limit, and they are difficult to obtain analytically, 
since T\ and Ti lie outside the Fermi liquid regime where 
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FIG. 4: Dependence of the temperatures T\ (filled symbols, 
solid lines) and T2 (open symbols, dashed lines) at which S(T) 
changes sign as a function of v 9 > 0.25 for U/T = 3, 6, 8 (left 
panel: quantum dot, right panel: magnetic impurity). The 
critical gate voltage, v 9 , beyond which no sign change occurs 
in S(T) at finite v 9 is indicated in the legend for each case. 



U/T 


T K (v g = o)/r 


Ti(v s = 0.25)/r 


T 2 (v 9 = 0.25)/r 


8 


2.64 x 10" 3 


0.044 (0.04) 


3.04 (1.39) 


6 


1.10 x 10~ 2 


0.056 (0.05) 


2.12 (1.00) 


3 


8.20 x 10~ 2 


0.13 (0.12) 


0.71 (0.41) 



TABLE I: Kondo temperature Tk (in units of T) at mid- 
valley v s = (symmetric point) for U/T = 8, 6, 3. Also shown 
are the temperatures Ti and Tb (in units of T) at which the 
thermopower, S(T), changes sign at v g — 0.25 (the smallest 
v 9 studied). The numbers in brackets are the corresponding 
temperatures for magnetic impurities. 



analytic calculations are possible. Estimates of these val- 
ues at the smallest gate voltage are close to the limiting 
values. They are tabulated in Table HI together with the 
relevant Kondo scales at mid-valley. Whereas the limit- 
ing values of T\ are comparable for both quantum dots 
and magnetic impurities, the limiting values of T2 for 
quantum dots are approximately twice larger than for 
magnetic impurities. By carrying out additional calcula- 
tions, using a finer grid of gate voltages, we determined 
the critical gate voltages v c g {U /T) : beyond which no sign 
change occurs (indicated in Fig. |3^-b) . For each value of 
U /r, we find that v g corresponds to entering the mixed 
valence regime, i.e. v g corresponds to a local level posi- 
tion Sd ~ — r/2 in the single channel Anderson model. 



C. Thermal conductance: K B (T) 

The electronic contribution to the thermal conduc- 
tance of a strongly correlated quantum dot, shown in 
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Fig. |3g-i, also exhibits interesting behavior: a cross- 
ing point at T w Y is found in the Kondo regime 
and for gate voltages approaching the mixed valence 
regime (Fig. [3g and inset). Such (approximate) cross- 
ing points are typical signatures of strong correlations 
and are well known in other contexts, including 3 He 
and heavy fermionsaa, dissipative two-level systems^ and 
doped Mott-insulators^. On entering the mixed valence 
and empty orbital regimes (Fig. [3ji-i) , two-peaks develop 
on either side of the crossing point (the lower peak being 
at T < 0.5Y and the upper one at T > 2Y). These quali- 
tative features in K e (T) can be related to A(ui,T), as in 
the case of S(T) (see also See. EE}- 



D. Moderate to weak correlations 

The effect of reducing correlations to a moderate value, 
U/Y = 3, is shown in Fig. [TT]of Appendix |B| the trends 
arc similar to those described above, with a significantly 
diminished Kondo regime. In particular, the evolution 
with gate voltage of G(T) is similar to that in the strongly 
correlated case (see Fig. [TTk-c) and the thermopower ex- 
hibits two sign changes as a function of temperature in 
the Kondo regime (Fig. [TTH) , with a rapid evolution to a 
single negative peak in the mixed valence and empty or- 
bital regimes (Fig. lllb -f). However, the crossing point in 
K c in the Kondo regime becomes less evident for moder- 
ate correlations (Fig.QTg), and, the two-peaked structure 
for K e in the mixed valence and empty orbital regimes is 
replaced by a single peak with a shoulder (Fig. Illh-i). 

These general trends, for correlated quantum dots, 
contrast with those for weakly correlated quantum dots, 
shown in Fig. [I2h-c of Appendix [B] for U/Y = 1 . These 
exhibit no sign change in the thermopower for any gate 
voltage v g > 0. Similarly, the thermal conductance for 
weakly correlated quantum dots shows no crossing point, 
exhibiting only a single finite temperature peak. 



E. High temperature asymptotics 

The FDM approach allows us to easily investigate the 
high temperature asymptotics of transport properties. 
As we discuss also in the context of dilute magnetic im- 
purities in Sec. IVII below, earlier transport calculations^ 
could not discern the highest temperature peak in S(T) 
(occurring at T « 6r for U/Y = 8, see Fig. [3Ji), nor 
the peak in the thermal conductivity (see discussion in 
Sec. IVII below). Here, we are able to do so. In addition, 
the numerical calculations recover the high temperature 
asymptotics of the transport properties: G{T) ~ 1/T, 
S{T) ~ 1/T and K C (T) ~ 1/T 2 for T > Y. Note that, 
for the Anderson model, the logarithmic corrections in 
the Kondo regime occur at intermediate temperatures 
T K « T « Y: the corrections at T > Y go over to 
the above power laws. 



VI. COMPARISON WITH DILUTE MAGNETIC 
IMPURITIES 

It is interesting to quantify the differences in the trans- 
port properties of quantum dots given by (OH with the 
analogous transport properties of dilute magnetic im- 
purities given by (|lH!13p . This is shown in Fig. [5] for 
the temperature dependence of transport properties in 
the Kondo, mixed valence and empty orbital regimes for 
U/Y = 6. 

A. Comparison of G(T) and p(T) 

In the Kondo regime, and for temperatures T << Y, 
the conductance of a quantum dot is a universal func- 
tion of T/T K , i.e. G(T)/G(0) = f(T/T K ) (e.g. see 
Ref. [47h . The same holds for the analogous quantity 
for dilute magnetic impurities, namely the resistivity, i.e. 
p(T)/p(0) = f{T/T K ) (e.g. see Ref. EH). Since G and 
p are different physical quantities, the functions / and 
/' are different and they cannot be made to coincide by 
using a common Kondo scale Tk (e.g. the Kondo scale 
defined in Eq. [T9|) . This is seen in Fig. [5^, which shows 
that the conductance curves for quantum dots are shifted 
in temperature, on a logarithmic scale, relative to the re- 
sistivity curves of magnetic impurities. The two functions 
/ and /' are rigorously identical only in the Fermi liquid 
regime T <C Tk- Experimentally, however, the accessible 
range of temperatures is that around T « Tk, say one 
decade below and one decade above Tk ■ For this region of 
temperatures, the two functions / and /' can be made to 
coincide by redefining them as new functions / and /', re- 
spectively, with different respective Kondo scales, T^ and 
T K such that f(T/T^ = 1) = f'(T/T K = 1) = 1/2, see 
Ref. |47|. In the mixed valence and empty orbital regimes, 
Fig. [SJd-c shows that the conductance of a quantum dot 
differs significantly from the resistivity of magnetic im- 
purities (with significant deviations at T > O.ir). In 
particular, the aforementioned finite temperature peak 
in the conductance of a quantum dot is absent in the re- 
sistivity of magnetic impurities. A signature of this peak 
in p{T) is seen at most in the Kondo regime at tempera- 
tures of order Y (see Fig. and is absent in the mixed 
valence and empty orbital regimes. These differences to 
the quantum dot case, arise, as described in Sec. IIIII 
due to the different way in which the spectral function 
appears in the respective transport integrals. These dif- 
ferences reflect also the absence of universality outside 
the Kondo regime. 

B. Comparison of thermopowers: S(T) 

In Fig. [5ji-f we see that, up to an overall sign change, 
due to A(u,T) appearing differently in the transport in- 
tegrals as explained in Sec. IIIII the thermopower of mag- 
netic impurities behaves in a qualitatively similar way to 
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that of a quantum dot, with two sign changes at Ti(v fl ) 
and T 2 (v g ) (shown in Fig. 0t>) and three extrema. In 
the Kondo regime, the position, T p , of the Kondo en- 
hanced peak in the thermopower of magnetic impurities 
is found to scale with Tx, just as for the quantum dot 
case (see Sec. lVB|) . A significant difference between S(T) 
for magnetic impurities and quantum dots is the much 
larger high temperature peak (at T > T2) for the for- 
mer in the Kondo regime (by as much as a factor 5, see 
Fig. [5ji). This difference holds to some extent also in 
the mixed valence regime (Fig. [5^). In the empty orbital 
regimes the thermopowers show a single peak at T sa T 
with a similar magnitude for both cases (Fig. |5j) . 

C. Comparison of K C (T) and K e (T) 

The electronic contribution to the thermal conduc- 
tivity of magnetic impurities, k c , shows significant dif- 
ferences to the corresponding thermal conductance of 
quantum dots, K e , sec Fig. [5g-i. For example, whereas 
K e exhibits interesting structure with either one (in 
the Kondo regime) or two (in the mixed valence and 
empty orbital regimes) peaks around T « T, n c only 
exhibits a single peak in all regimes, and this peak 
occurs at a much larger temperature T ^> T. The 
reason for the latter difference is the following: the 
main contribution to the thermal conductance and ther- 
mal conductivity come from the integrals I2 and M 2 in 
Eq. ([7]) and Eq. (|14[) .respectively, which involve integrals 
of -(df/du)uj 2 A(u],T) and - (d f / duj)uj 2 / A(uj , T), re- 
spectively. For the former, the peaks in A(u>, T) at sa and 
Ed + U result in a peak in the integrand at |w| > U/2 > T, 
whereas in the latter, the dips in l/A(u,T) at Ed and 
Ed + U shift the main contribution to the integral to much 
higher energies |w| 3> U/2. Correspondingly, the temper- 
ature of the peaks in K e in the former are at T < U/2 
and for the latter are atT^J7/2>r, in agreement 
with the numerical results. The existence of two peaks 
in K e in the mixed valence and empty orbital regimes as 
opposed to a single peak in the Kondo regime is also eas- 
ily explained: the two peaks reflect the sampling of the 
two incoherent features at Ed and Ed + U in u 2 A(ui,T) ap- 
pearing in the moment I 2 for K e . In the Kondo regime, 
these excitations lie close to each other and only one peak 
results. Similarly, the single peak in K e for all regimes re- 
sults from the strong suppression of the above incoherent 
excitations in 1/A(uj,T) appearing in the moment M 2 . 



VII. FIGURE OF MERIT, POWER FACTOR 
AND LORENZ NUMBER 

A measure of the thermoelectric efficiency of a quan- 
tum dot device is the dimcnsionless figure of merit de- 
fined by ZT = GS 2 T/(K C + K ph ), where K ph is the 
phonon contribution to the thermal conductance. Hence 
for high efficiency, one requires either large S or small 



total thermal conductance relative to electrical conduc- 
tance or both conditions simultaneously. A calculation 
of ZT for quantum dot systems would therefore require 
knowledge of the material specific phonon contribution 
to the thermal conductance K^. Similarly, for magnetic 
impurity systems a calculation of the dimcnsionless figure 
of merit, ZT = aS 2 T /(k c + Kph), would require knowl- 
edge of the material specific phonon contribution to the 
thermal conductivity K p h- This is outside the scope of 
the present paper, so instead we show in Fig. |6K-c re- 
sults at U/T = 8 for the quantity ZT = GS 2 T/K C , for 
quantum dots, and ZTq = aS 2 T j K e , for magnetic im- 
purities (with the latter being depicted on the negative 
axis for clarity). In addition, we also show in Fig. [6jl-f 
an appropriate rescaled power factor (PFq = S 2 G/Gq for 
quantum dots, and PFq = S 2 o /oq for magnetic impurity 
systems). This is another useful measure of an efficient 
thermoelectric system, by-passing lack of knowledge of 
the total thermal conductances (conductivities). 



A. Figure of merit 

Since, in the low temperature limit, the thermopower 
of the Anderson model vanishes linearly with tempera- 
ture in all regimes, a significant figure of merit is found 
only at finite temperature, as seen in Fig. [6K-c. 

In the Kondo regime, Fig.[6K, enhanced regions of ZTq 
are found in four temperature regions, (i), at T ~ T p , 
due the low temperature Kondo enhancement of the ther- 
mopower, however the magnitude of ZTq is tiny (see inset 
to Fig. [6^), (ii), at temperatures of order T in the region 
Ti <T < T 2 , where ZT can be of order 0.2-0.3 for both 
quantum dots and magnetic impurities, (iii), at temper- 
atures T > T2, with enhancements comparable to those 
for region (ii), and, (iv), in the asymptotic region T»r, 
where ZTq saturates to a finite value which is larger for 
quantum dots than for magnetic impurities (discussed 
below) . 

The behavior of the figure of merit in the mixed valence 
and empty orbital regimes is complicated, see Fig. [BJd-c. 
In the mixed valence regime, significant enhancements 
are found, for quantum dots, at temperatures some- 
what below r, see Fig. [6Jd, and in the asymptotic regime 
T 3> r. Similar enhancements arc found also for the 
empty orbital case (Fig. [6]:). For magnetic impurities, 
similar enhancements to quantum dots are found on tem- 
perature scale of order T, but at T 3> T the enhancements 
are much smaller than for quantum dots (see Fig. |SJ> 
c). The latter effect is due to the much larger thermal 
conductivities (even at higher temperatures) of magnetic 
impurities as compared to those of quantum dots (see 
discussion above and Fig.[Sg-i). 
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Kondo Mixed valence Empty orbital 




T/T T/T T/T 

FIG. 5: Comparison of transport properties for quantum dots (solid lines) and magnetic impurities (dashed lines) in the strongly 
correlated regime U/T = 6. (a-c): normalized electrical conductance G(T)/Go (quantum dot) and normalized resistivities p/po 
(impurity), where po is defined in (|16p . (d-f): thermopower S(T) for quantum dots and impurities (inset to (d) shows the 
low temperature Kondo peak in the thermopower in more detail), (g-i): electronic contribution to the thermal conductance 
ciKe/Go (quantum dot) and thermal conductivity C2K e po (impurity), rescaled by Go and po respectively, so that the same 
unit (K) applies to both cases. The numerical factors C\ = 10 2 and C2 = 5 x 10~ 4 are included for clarity of presentation. 
The inset in (g) for K e (T)/Go shows the crossing point in K e (T)/Go at T/T ~ 0.6 in more detail and the evolution of the 
second peak in the thermal conductance of the quantum dot as the mixed valence regime is approached (with inclusion of 4 
additional gate voltages v 9 = 1.50, . . . , 2.25 in the Kondo regime). The range of gate voltages is v 9 = 0.25, . . . , 1.25 (Kondo), 
v g = 2.25, . . . , 3.75 (mixed valence), and v 9 = 4.0, . . . , 5.75 (empty orbital). 



B. Power factor 

The power factor PFq is enhanced in the same regimes 
(i)-(iii) as the figure of merit, see Fig. [6^-f, but vanishes 
as 1/T 3 in the limit T » T (using the asymptotic be- 
havior of S and G from Sec. IV E[) . In the Kondo regime, 
PFq exhibits a much larger peak above T% for magnetic 
impurities as compared to quantum dots (see Fig. [SJl). 
This reflects the observation made above (Sec IVI Bp that 
the highest temperature peak in S(T) at T > Ti for mag- 
netic impurities is significantly enhanced as compared to 
that of quantum dots. For quantum dots, the main en- 
hancement in PFq in the Kondo regime is in the range 
T < T<z. These trends differ little from those observed 
in the mixed valence and empty orbital regimes for both 
quantum dots and magnetic impurities (see Fig. |B]>f). 



C. Wiedemann- Franz law and Lorenz number 

We comment on the enhancement of ZTq in the re- 
gion T > T. which can result in ZTq > 1 (e.g. in the 
mixed valence and empty orbital cases). This enhance- 
ment reflects a violation of the Wiedemann- Franz law at 



T r. The latter states that the thermal conductance 
(conductivity) is proportional to the electrical conduc- 
tance (conductivity) multiplied by temperature, i.e. that 
the Lorenz number L(T), defined for quantum dots by 

L(T) = K C (T)/TG(T), (21) 

and for magnetic impurities by 

L(T) = n e (T)/Ta(T), (22) 

is independent of temperature and takes on the universal 
value L = n 2 k 2 B /3e 2 . Since, ZT = S 2 /L{T), a sig- 
nificant reduction of L(T)/Lq can result in an enhance- 
ment of ZTq. In Fig. |5g-i we see that L(T)/Lq is much 
suppressed at T ^> T, thereby allowing for significant 
enhancements in ZTq in this limit. This enhancement 
is seen for all regimes, especially for the mixed valence 
and empty orbital regimes. We note, however, that from 
Fig.[6g-i (and inset), the Wiedemann- Franz law is, on the 
whole, reasonably well satisfied at temperatures T«T, 
and becomes exact in the Fermi liquid rcgimei^(for other 
violations of the Wiedemann- Franz law see Ref |58[). 
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FIG. 6: Temperature dependence of, (a-c), the "figure of merit", ZTq, (d-f), the "power factor", PFq, and, (g-i), the Lorenz 
number ratio L{T)/Lq for a quantum dot (solid lines, positive y-axis) and for magnetic impurities (dashed lines, shown on 
the negative y-axis for clarity). Results are for U/T = 8 and the same range of dimensionless gate voltages as in Fig. [3] For 
quantum dots we define ZT = GS 2 T/K C , PF = S 2 G/G and L(T)/L = K B /GT with L = 7r 2 fc|/3e 2 . The corresponding 
quantities for magnetic impurities are defined by ZTq = uS 2 T/k c , PFq = S 2 a/ao and L(T)/Lo = (n e /aT) where a = 1/p and 
cto = 1/po • Arrows indicate the evolution of the transport quantities with increasing v 9 > 0. Insets in (a) and (d) for the 
Kondo regime, show the low temperature peak in ZTq and PFo in the vicinity of Tk(v 9 ) > 2.6 x 10 r (Table [I]). The inset 
in (g) for the Lorenz number shows the deviations from the Wiedemann-Franz law in the region around Tx(v 9 ). In this inset, 
L(T)/Lq is shown on the positive y-axis for both impurity and quantum dot cases. 



VIII. UNIVERSAL SCALING FUNCTIONS FOR 
THERMAL TRANSPORT THROUGH 
QUANTUM DOTS 

By analogy to the scaling properties of the electrical 
conductance G(T)/G(0) = f(T/T K ), where / is a uni- 
versal function of t = T/Tk in the Kondo regime, it 
is interesting to establish to what extent such scaling is 
present in the thermopower, S(T), and the thermal con- 
ductance, K e (T), of strongly correlated quantum dots. 
We investigate this here, for U/T = 8 and for values of 
the gate voltage in the Kondo regime (see Fig. [7]). 

In the Fermi liquid regime, T <C Tk, we havei^ 



S(T) 



TTjT 



cot(7rnd/2). 



(23) 



Scaling can therefore be expected for S(T)/T, once 
the occupancy (and gate voltage) dependent factor 
cot(7rnd/2) is scaled out. In the above, 7 ~ 1/Tk 
is a measure of the inverse Kondo scale and can be 
extracted^ from the numerical value of limr^o \S(T)\/T 
using Eq. ([2"3"T) and the calculated values of n d from 
Fig. [Q We sec from Fig. that s(T/T K ) = 
\e\S(T)/ir^T cot(irnd/2) does indeed scale with 7T ~ 



T/Tk for a range of gate voltages in the Kondo regime. 
This scaling extends up to temperatures comparable to 
Tk with significant deviations setting in above this tem- 
perature scale. This is not surprising given the fact 
that the thermopower is a highly sensitive probe of the 
particle-hole asymmetry in the spectral density. 

For the thermal conductance, K e (T), we expect from 
the Wiedemann-Franz law, K c /T ~ G(T), to see a scal- 
ing in K C (T)/T similar to that in G(T)/G(0). This is 
confirmed in Fig.Eb which shows K C (T) juT versus T jT K 
for several gate voltages in the Kondo regime, where a is 
defined by 



lim 

T-J-O 



T 



and T K is a Kondo scale defined by 



K C (T = T K ) 



T 



K 



n 
2' 



(24) 



(25) 



We see that, for U/T = 8, K c (T)/aT = g{T/T K ) is a 
universal function of T /T K for temperatures extending 
up to at least lOOTg., just as G(T)/G(0) = f(T/T K ) is 
a universal function of T/Tk for temperatures extending 
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FIG. 7: (a) S(T)/T for a quantum dot, scaled by its limit- 
ing low temperature absolute value 717 cot(7rnd/2) /|e| versus 
7T ~ T/Tk for (7/r = 8 and a range of gate voltages in 
the Kondo regime. Here, 7 ~ 1/Tk^. (b) the thermal con- 
ductance K e (T), scaled by aT, versus T/T^, for parameters 
as in (a), with a defined in (|24[) and a Kondo scale de- 
fined in (|25[) . The electrical conductance G(T)/G(0) versus 
T/T K at v g — 0.25 is also shown. The Kondo scale defined 
by G(T = Tk) = G(0)/2 is related to T K by T K w 1.9T&. 

up to at least lOOT^. Increasing U/T, and thereby re- 
ducing Tk /r allows universality to extend to still higher 
temperatures. Suppressing charge fluctuations, e.g. by 
working within a Kondo model, allows these universal 
scaling functions to be defined for all temperatures. Note 
also, that although these universal functions / and g have 
a similar functional dependence on T/Tk and T/T K re- 
spectively, they are shifted relative to one another on 
an absolute temperature scale. The difference between 
g and / for temperatures around T = T K w 2\-/1.9 
accounts for the violation in the Wiedemann-Franz law 
on this scale, as noted previously (see inset to Fig. |Bg). 
The Wiedemann-Franz law is only satisfied exactly in 
the Fermi liquid regime T -C Tk ■ One can collapse G(T) 
onto K e (T)/T by scaling the temperature axis of the for- 
mer by T k /Tk- In the universal regime T « T, small 
deviations between g(t) and f(t) arise for t < 1 and t > 1. 

For dilute magnetic impurities, our conclusions for 
scaling in the Kondo regime are essentially the same as 
those above for quantum dots (see also Ref. Il5[ ). 

IX. GATE VOLTAGE DEPENDENCE OF 
TRANSPORT PROPERTIES 

The gate voltage dependence of transport through 
quantum dots and magnetic impurities is shown for sev- 
eral representative temperatures in Fig. [5J for the elec- 
trical and thermal transport and in Fig. [9] for the ther- 
mopower. For magnetic impurities, different v g should 
be understood as corresponding to changes in the local 
level position Sd relative to the Fermi level, as invoked 



by the application of pressure (either chemical via dop- 
ing or hydrostatic). In some rare earth systems^ the 
application of pressure has been shown to tune the mag- 
netic impurities from the Kondo to the mixed valence and 
empty orbital regimes. Throughout this section we de- 
note by T\ and T2 the minimum temperatures of Ti(v s ) 
and T 2 (v g ) in the limit v g ->• (see Figg]and Tabic [§. 
We show that the different behavior in the gate voltage 
dependence of the thcrmopower, at fixed temperature T, 
can be classified in terms of the relative value of T to T\ 
and T2. 



A. Gate voltage dependence of G and K c 

The gate voltage dependence of the electrical and ther- 
mal conductance of quantum dots is shown in Fig. [5^,-c. 
The former exhibits, for T ^> Tk, Coulomb blockade 
peaks at Ed = ±U/2T with a suppression of the con- 
ductance in the mid- valley region around v g = 0. On 
decreasing the temperature, the Kondo effect becomes 
operative resulting in an enhancement of the conduc- 
tance in the region between the Coulomb blockade peaks 
(see Fig. |8]i). This picture is well known. At T > T K , 
the thermal conductance of quantum dots also exhibits 
Coulomb blockade peak o 60 ' 61 , not directly evident in the 
plot of K e versus gate voltage (Fig. [5b), where only weak 
signatures of these are discernible. They become clearer 
in the gate voltage dependence of K e /T which by the 
Wiedemann-Franz law (see Sec. lVII C|) is proportional to 
the electrical conductance G(T), as seen in Fig. [SJ;. Dif- 
ferences between Fig. [8^, and Fig. [5J: indicate the degree 
of deviation from the Wiedemann-Franz law. These de- 
viations are largest at T T for all gate voltages, as 
previously observed in Fig. Hg-i. 

The same observations, using a different terminol- 
ogy, can be made for the case of magnetic impurities 
in Fig. [SU-f: the valence fluctuation peaks at T ^> T 
arc now seen in the resistivity of mixed valence impu- 
rities, whereas Kondo impurities (v g « 0) have a small 
resistivity at T > T and a large unitary resistivity at 
T = 0. The behavior of the thermal conductivity k is 
similarly understood in terms of the Wiedemann-Franz 
law k c (T) ~ aT with a = 1/p, as seen by comparing 
Fig. [11 and Fig. [SJ. 



B. Gate voltage dependence of S 

The gate voltage dependence of the thcrmopower of 
quantum dots is shown in Fig.|9^i. The particle- hole sym- 
metry about v g = (see Sec. HJ) implies S- v = —S v 
at all temperatures, so we only discuss v g > 0. We focus 
mainly on the Kondo regime, v g < Vg, and discuss the 
remaining gate voltages by reference to Fig. [3] There are 
three main types of behavior, characterized by the follow- 
ing temperatures, T, relative to T\ and T2: (i), T < Ti, 
as exemplified by T = O.Oir, (ii), T x < T < T 2 as ex- 
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FIG. 8: Gate voltage dependence of electrical and thermal transport properties of quantum dots (left panels (a-c)) and magnetic 
impurities (right panels (d-f)) for U/T = 8 and at several representative temperatures T/T (labeling the individual curves in the 
figures). Vertical dashed lines in the top panels delineate the Kondo regime from mixed valence and empty orbital regimes, (a): 
The normalized conductance G(T)/Gq ( (d): the normalized resistivity p(T)/po for the impurity case), and, (b), the rescaled 
thermal conductance K c (T)/Go ( (e): the rescaled thermal conductivity K c (T)po for the impurity case). The bottom panel, 
(c), shows the dimensionless ratio (K c /T)/(LoGo) for the quantum dot (panel (f) shows the analogous dimensionless ratio for 
magnetic impurities (k c /T) _1 /{pq/Lq)). Deviations from G(T)/Go in (a) reflects deviations from the Wiedemann-Franz law 
(similarly, deviations of (K e /T) _1 /(po/I/o) in (f) from p(T)/po in (d) indicate deviations from the Wiedemann-Franz law for 
the impurity case). For clarity not all temperatures are shown in (b) and (e). 
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FIG. 9: Gate voltage dependence of the thermopower for the 
strong correlation case U/T = 8 at several characteristic tem- 
peratures T/T for, (a), the quantum dot case, and, (b), the 
impurity case. The vertical dashed lines delineate the Kondo 
regime (around v 9 = 0) from the surrounding mixed valence 
and empty orbital regimes. 



emplified by T = O.ir and T = T, and, (iii), T > T 2 , 
as exemplified by T = 10r. In case (i), the Kondo 
resonance is asymmetric about the Fermi levels, lying 
slightly above it for v g > 0. The slope of the spec- 
tral density at lu = is positive, resulting by Eq. ([20)1 



in a negative thermopower, as observed for T = 0.0 IT. 
The same holds, at still lower temperatures, T <C Tjc, 
where Fermi liquid theory^ gives the explicit expression 
(1231 . Case (ii), T\ < T < T 2 , is the most interesting for 
quantum dots, for several reasons: first, this temperature 
range is experimentally accessible since for U/T = 8, we 
have Ti = 0.044r and T 2 = 3.04r. Second, there is an 
overall sign change in S(v g ), relative to case (i), for a fi- 
nite range of gate voltages (see Fig.[S^). Third, a further 
sign change occurs at finite v g > 0, and, fourth, the ther- 
mopower is large enough for a significant range of gate 
voltages to enable its measurement. The sign change at 
a finite gate voltage occurs when S(T\ < T < T 2 ) as 
a function of v g in Fig. [3Ji reaches the value zero and 
becomes negative. 

For gate voltages outside the Kondo regime, the ther- 
mopower, as a function of gate voltage, S(v g ), either ap- 
proaches zero at v 9 > 1 (as happens for T — 0.1 in 
Fig-Hk) or does not saturate for v g ^> 1 (e.g. for T = T in 
Fig. |Hfi). In terms of Fig. [3f (see the arrows), the former 
occurs for temperatures to the left of the minimum in S 
in Fig.[5f, and the latter occurs for the opposite case. The 
latter case is half-way to case (iii), T > T 2 ^> T, which 
exhibits a thermopower approximately linear in gate volt- 
age, with no sign change at any v g > 0. This is similar 
to the "sawtooth" behavior of S(v g ) found for multi- 
level quantum dots weakly coupled to leads at T ^> T 
in Ref. l62ll63l (for related experimental work see Ref. l64h - 
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l69h . The behavior of the thermopower in multi- level open 
quantum dots has also been investigated^— . 

The same classification (i)-(iii), as for quantum dots, 
can be used to explain the local level dependence of the 
thermopower of magnetic impurities shown in Fig. ^jp. 



X. CONCLUSIONS AND DISCUSSION 

In this paper we investigated the thermoelectric prop- 
erties of strongly correlated quantum dots, described by 
the single level Anderson impurity model connected to 
two conduction electron leads. For this purpose, we used 
Wilson's NRG method and calculated the local Green's 
function and transport properties by using the full den- 
sity matrix approach^. Since this approach builds into 
the density matrix all excitations obtained in the NRG 
approach, it is particularly well suited to finite temper- 
ature transport calculations, allowing us, for example, 
to investigate also the high temperature asymptotics of 
transport properties. 

For strong correlations and in the Kondo regime, the 
thermopower exhibits two sign changes, at tempera- 
tures Ti(v s ) and T 2 (v 9 ) with T x < T 2 . Wc found that 
T\ > T p (v g ) « Tfs-(vg), where T p (v g ) is the position of 
the Kondo induced peak in the thermopower, r^(v g ) is 
the Kondo scale, and T 2 = 0{T). The loci of Ti(v s ) and 
T2(v g ) merge at a critical gate voltage v 9 = Vg(U /T), be- 
yond which no sign change occurs. We determined v g for 
different U/T finding that Vg coincides, in each case, with 
entry into the mixed valence regime. No sign change is 
found outside the Kondo regime or for weak correlations, 
U/T < 1. Thus, a sign change in S(T) at finite v g is a 
particularly sensitive signature of strong correlations and 
Kondo physics. This effect could be measurable in quan- 
tum dots, as it manifests itself in an overall sign change 
in S for a finite range of gate voltages on increasing tem- 
perature T from below T\ (v 9 — > 0) to values in the range 
T\ < T < T 2 = 0(T), which is an accessible range since 
Ti > T K . 

The results for quantum dots were compared also to 
the relevant transport coefficients of dilute magnetic im- 
purities in non-magnetic metals: the electronic contribu- 
tion, k c , to the thermal conductivity, the thermopower, 
S, and the impurity contribution to the electrical resis- 
tivity, p. As regards the temperature dependence of the 
respective transport quantities, we find, in the mixed va- 
lence and empty orbital regimes, two peaks in K e (T) as 
compared to a single peak in k c (T). Similarly, G(T) ex- 
hibits a finite temperature peak on entering the mixed va- 
lence regime, whereas such a pronounced peak is absent 
in p(T), even far into the empty orbital regime. As for 
quantum dots, we find that the low temperature Kondo 
peak position in the thermopower of magnetic impurities 
scales with Tk- We compared and contrasted the figure 
of merit, power factor and the extent of violation of the 
Wiedemann-Franz law in quantum dots and dilute mag- 
netic impurities, finding enhanced figures of merit at tem- 



peratures where the Wiedemann-Franz law is strongly 
violated. Finally, we clarified the extent of scaling, as a 
function of T/Tk, in the thermopower and thermal con- 
ductance of quantum dots in the Kondo regime. 

We comment on a recent experiment in Ref. [2C3 which 
we believe shows evidence of Kondo correlations in the 
thermopower of a strongly correlated quantum dot. In 
this experiment, the thermo voltage across a Kondo cor- 
related quantum dot is investigated as a function of gate 
voltage and lattice temperature. This can be compared 
to our S(v g ) in Fig. [9^. The gate voltage Ve in Ref. 
is related to our dimcnsionless gate voltage, v s , via 
— eVg ~ £d ~ v g , i.e. Ve ~ — v ff . Mirror reflecting our 
results for S(v g ) in Fig.^ about v g = allows a qualita- 
tive comparison with the experimental measurements in 
Ref. H3. Using the experimental estimate T = 0.35 meV 
from Ref. [2p| . we can translate the four experimental tem- 
peratures T L = 0.07K, 0.25 K, LOOK and T L = 1.46K 
at which the thermopower was measured into our the- 
oretical temperatures in units of T. We assume strong 
Coulomb correlations on the dot U/T = 8 and show 
the results for the thermopower in Fig. [10J From Ta- 
ble |H the lowest experimental temperature corresponds 
to T <C Ti, the next lowest temperature (T = LOOK) lies 
close to T = Ti , where the thermopower changes sign in 
the Kondo regime, and the highest two temperatures lie 
between Ti and T%. The lowest temperature measured, 
T = 0.07 K, indeed shows a positive thermopower above 
mid- valley, in agreement with our results for T < T\. 
Upon increasing the temperature, the experiment shows 
a sign change of the thermovoltage for a finite range of 
gate voltages (relative to mid- valley) , which is consistent 
with our prediction of such a sign change in the Kondo 
regime for T 2 > T > Ti. The onset, with increasing 
temperature, of an additional oscillation in S(y g ) about 
v g = in the experiments is therefore consistent with 
our results. The experimental data deviates from our 
calculated thermopower in the mixed valence and empty 
orbital range of gate voltages v g 3> 1, with the theoret- 
ical results showing a much larger thermopower in this 
region of gate voltages. These deviations are expected at 
v 9 > 1 since additional levels present in real quantum 
dots, but absent in our model, start being populated. 
This significantly influences transport through the quan- 
tum dot. Qualitatively, however, we are able to interpret 
these experiments on the thermopower of Kondo corre- 
lated quantum dots for gate voltages v g ks 0. For a more 
quantitative comparison to theory, further investigations 
are needed. 

Calculations for single quantum dots and dilute mag- 
netic impurities, are a starting point for dealing with a 
finite density of quantum dots, such as self-assembled 
quantum dots, or for a finite concentration of magnetic 
or mixed valence impurities in bulk (e.g. for TI impu- 
rities in PbTe in Ref. @). The transport properties of 
such systems, modeled by a random distribution of An- 
derson impurities, will be determined by (|11|13|) subject 
to the charge neutrality condition niUd + n c = n, where 
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rid is the occupancy of the dot (impurity) , n c is the occu- 
pancy of the relevant conduction band and n is the total 
electron filling. Coupled with material specific electronic 
structure information and the effects of phonons, such 
calculations, will be important for understanding the po- 
tential of materials such as self-assembled quantum dots 
or PbTei-^TL; systems for thermoelectric applications. 
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Appendix A: Reduction of two-channel Anderson 
model to a single-channel Anderson model 

The reduction of the single-level two-lead Anderson 
model ([1} for a quantum dot, to a single-channel model 
is, in general, approximate, but as we show here, the ap- 
proximation is very good (or even exact). One notices 
first, that the d-statc of the quantum dot in ([TJ only cou- 
ples to the even combination t L c Lka + ^RC-Rka of the lead 
electron states. By using the following canonical trans- 
formation 



(o) lead states, as follows 

H = ^ e eka a\ kiy a eka + ^ £d d ° d,J + Un dt n di 

ka a 

+ t^(al ka d a + h.c.) + H + H pot . (A3) 



ka 

Here, e eka = 

Sfcer e okaCl oka .a krT 



(£Lkcrt L + £ Rkatft) / 't 2 , H Q = 

is the Hamiltonian for the odd 
lead electrons with e oka = (e Lka t 2 R + e Rka t 2 L )/t 2 , 
and H poi = Y, ka - U k (a\ klT a oka + h.c.) is a potential 
scattering term between even and odd lead electrons. 
Hence, the odd lead electrons do not couple to the 
dot directly, but only indirectly via the potential 
scattering term. The magnitude of this is given by 
U k ° = {shko - £Rka)tLtR./t 2 , which is vanishingly small 
at low energies. Moreover, it vanishes identically for 
degenerate leads e L ka = £Rka- The calculations we 
report in this work, using the single channel Anderson 
model 

H = ^ e eka al k(J a ekc , + ^ e d d\d a + U n d ^n d ^ 

ka a 

+ tJ2(4ka^ + h-c), (A4) 



.. Since t 2 = t 2 L + t 2 R , 
TiNpt 2 of the single chan- 



are therefore a very good approximation, even in general 
to those obtained from the two-lead model ([T]) and identi- 
cal to those for the case e i, ka = e ru 
the hybridization strength V 
nel model is seen to be the relevant single-particle broad- 
ening, Tl + Fr, of the two- lead model ([T]). In this paper 
we follow the convention in the quantum dot commu- 
nity of using the full-width at half-maximum, T = 2T, 
as the unit of energy. Finally, we note, that a reduc- 
tion to a single-channel model is, in general, not possible 
for multi-level or double quantum dots attached to two 
leads^-^. 



Appendix B: Results for moderate and weak 
correlations 



Fig. [TT1 shows results for a moderately correlated quan- 
tum dot, U/T = 3 exhibiting the same trends as those 
found for the strongly correlated case U/T — 8 (Fig. [3]). 

For completeness, we show an example of transport 
through a weakly correlated quantum dot with U/T = 1 
in Fig. [12] In this case, the thermopower remains neg- 
ative for all gate voltages v g > (Fig. [T2b). Similarly, 
the thermal conductance exhibits only a single peak for 
all gate voltages (Fig. W7b). 



t a eka 

i ®>oka 



thChka + t R CR k(T 
thC Rka — tRCLka, 



(Al) 
(A2) 



noting that normalization of even/odd states implies t 2 



t R , we can rewrite (fl]) in terms of even (e) and odd 



Appendix C: Green's functions within the FDM 
approach 

In this appendix we give an alternative derivation of 
the finite temperature Green's function within the FDM 
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FIG. 11: Temperature dependence of, (a-c), the normalized 
electrical conductance G/Go, (d-f), the thermopower, S, and, 
(g-i), the normalized thermal conductance, K e /Go, multiplied 
by a factor 10 2 for clarity of presentation , as function of T/r, 
in the moderately correlated regime U/T — 3 and a range of 
dimensionless gate voltages, v 9 = (ed + U/2)/T > 0, in the 
Kondo (first column), mixed valence (second column) and 
empty orbital (third column) regimes. The range of v 9 is 
indicated in the top panels for each regime and the increment 
used was 0.25. Arrows indicate the evolution of the transport 
quant 
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FIG. 12: Temperature dependence of, (a), the normalized 
electrical conductance, G/Go, (b), the thermopower, S, and, 
(c), the normalized thermal conductance, K e /Go, multiplied 
by a factor 10 2 for clarity of presentation, as function of T/V, 
in the weakly correlated regime U/T = 1 and a range of di- 
mensionless gate voltages, v 9 = (ed + U/2) /T. The range of v 9 
is indicated in the top panel and the increment used was 0.25. 
Arrows indicate the evolution of the transport quantities with 
increasing v 9 . 



approach of Weichselbaum and von Delft22. A concise 



derivation, implementing arbitrary abelian symmetries, 
has also been given in Ref. HH. We consider a general 
fermionic retarded Green's function 

Gab® = -i6(t)([A(t),B)+) 

= -i6{t)Tr[p{A{t)B + BA(t))], 
where A, B are fermionic operators, e.g. for the d-level 
Green's function of our quantum dot A = d a and B = S a . 
The trace is evaluated for an appropriate density matrix p 
by using the complete set of states introduced by Anders 
and Schiller—. These consist of the set of states \lem) = 
|Zm)|e) obtained from the eliminated eigenstates, \lm), of 
H m , and the degrees of freedom, denoted collectively by 
e, of the sites i = m + 1, . . . , N, where N is the longest 
chain diagonalized. The retained low energy states of H m 
are denoted \km), and \kem) = |fcm)|e) extends these to 
the Hilbert space of Hn by the additional environment 
degrees of freedom e of the sites i = m + 1, . . . , N. The 
eigenstates (retained and eliminated), \p = (k,e)m), and 
eigenvalues, E™ = ^ k l y of H m satisfy H m \pm) = E™\pm). 
Completeness of the states \lem) is expressed by^i 



1 = \lem')(lem'\ 

m'-mo+l le 



(CI) 



where niQ is the last iteration for which all states are 
retained. For iterations m > m , the set of states \p = 
(k, l)m) consists of both retained (k) and eliminated (/) 
states. The following decomposition of (|C1[) is useful^: 



^2\lem'){lem'\ = l+ + 1" (C2) 

m 1 — mo + 1 le 



J2 J2\ lem '^ lem '\ 



(C3) 



m'— mo + 1 le 
N 



= \lem!) (lem'\ 

m' — m-f-1 le 

= \kem) {kem\ 



(C4) 



ke 



where the last equation follows from the fact that the 
Hilbert space of retained states at iteration m (sup- 
plemented by the degrees of freedom e for sites m' = 
m + 1, . . . , N) spans the same Hilbert space as all elimi- 
nated states from all subsequent iterations. By using the 
decomposition of unity (|C1[) twice within the trace in the 
expression for Gab®, the following Lehmann represen- 
tation can be found for this Green's function^ 
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G AB (t) 



AB — 



AB — 



'AB 



AB ' ^ifl b ,4Bi 
TV 

- i6> (*) E 2 [ et(Br ~ B!7)t (^m|A|/'e'm)(/'e'm|Bp|Ze?7i) + e^'^ (lem\A\l'e'm)(l'e'm\pB\lem) 

m=m + l le,l'e' 
N-l 

6(t) E [e t{Er ' E ' kn)t ( k e'm\Bp\lem)(lem\A\ke'm) + e t ^- E ^ )t {lem\pB\ke'm)(ke'm\A\lem) 



m=m + l le,ke 
N-l 

m E E [ e 

m=mo + l lem,ke 



l{E ™- Er)t {lem\Bp\ke'm)(ke'm\A\lem) + e l{EI "- E ^ )t {ke'm\pB\lem){lem\A\ke'm) 
I 



where the double sum over to, to' (coming from two appli- 
cations of (|Cip ) is decomposed into contributions to' = to 
(first term), to' > to (second term) and to' < m (third 
term). In the last two terms, use has also been made 
of (|C4|) . In the time evolution e lHt \pem),p = {k,l), we 
have made use of the NRG approximation H « H m , so 



\pem) 



HE" 



\pem). Peters et al^ evaluated 



that 

the above expression for the Green's function by using 
an approximate density matrix p^, defined by the elim- 
inated states of the longest chain diagonalized, i.e. 



1 



Pn 



N 



I 



\lN)e 



(IN\ 



(C5) 



where Pn = 1/A;bTjv (or Tjy) is chosen appropriately^ 
to ensure that Zn{Tn) is a good approximation to the 
partition function of the infinite system at tempera- 
ture T = Tm- Note also, that since N is the last it- 
eration, all states in the above expression arc consid- 
ered as eliminated states in order that (|Cl[) be satisfied. 
This procedure can be repeated for each chain length 
to = N, N — 1, . . . , toq + l , using a density matrix 



1 



i(T m ) ; 



\lm)e 



-f> m E? 



(C6) 



to obtain shell Green's functions G m (w),m = N,N — 
1, . . . , too + 1 defined at a corresponding set of temper- 
atures T m (or fi m = 1/ksTm) (for clarity we hence- 
forth omit the subscript AB for Gab)- Since these shell 
Green's functions, G m , contain only excitations of or- 
der the characteristic scale, io m , of H m , or larger, the 
Green's function G m (ui) can only be evaluated at fre- 
quencies lo > T m . Information at uj <C T m is not avail- 
able. This restriction is overcome by the FDM approach 
that we now describe. 

Weichselbaum and von Delft^ evaluated the above 
Green's function by using the FDM of the system made 
up of the complete set of eliminated states from all it- 
erations to = mo + 1,...,N. Specifically, the FDM is 
defined by 



N 



E E i^ em ) 



-0EV 



m=m + l le 



Z{T) 



■ (lem\ 



(C7) 



where Z(T) is the partition function made up from the 
complete spectrum, i.e. it contains all eliminated states 
from all H m , to = mo + 1, . . . , N. Consequently, evaluat- 
ing the Green's functions by using the above FDM, allows 
an arbitrary temperature T to be used for all frequencies 
w, and, in particular, allows accurate calculations to be 
carried out at lu <C T. 

Consider the following density matrix for the m'th shell 
(defined, however, in the Hilbert space of ifjy): 



le 



\lem) — = (lem\ 



(C8) 



Normalization, Tr[p m ] = 1, implies 



4 N- m Zrn (cg) 

Z m 



where Z m = J2i er m " '■ Then the FDM can be written 
as a sum of weighted density matrices for shells to = 
to + 1,..., TV 



iV 



E w mPm 
m— mo + 1 



iY 



iN-m 



(CIO) 
(Cll) 



m— mo + 1 



The calculation of the weights w m is outlined in Sec. IC 11 
Substituting p = J2 m ' w m'Pm' m t° the above Lehmann 
representation for G(t) and Fourier transforming yields 
G{u) = Em/ w m ,(& m ,{u) + G^,(u) + G^,{u)). The first 
term, G l m , , is easily evaluated by using the orthonormal- 
ity of the eliminated states (l'e'm'\lem) = Sui6 ee >5 mm i , 
orthonormality of environment degrees of freedom in 



(I 



em\A\l'e'm'} 



= 5ee<A™, , with Aft, = (lm'\A\l'm') and 



the trace over the N — to' environment degrees of freedom 



E 



(JV-r 
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to obtain (for m! — mo + 1, . . . , N) 



lu + Ef - Ef + id ' 



The second term, / , is also easily evaluated and results 
for m' = ?7io + l, . . . , N — 1 (the iV'th term vanishes, as all 
states are counted as eliminated states at this iteration) 



An ' ik 



lu + Ef - Ef + id 



-f)EV 



lu + Ef - Ef + id 



The third term, G l ™, , takes the form 

(kem\p m i B\lem) 



lo + ET 1 - E" 1 + id 

lek 1 K 



E A kl 



{lem\B p m i \kem) 
lu + E™ - E™ + i5 



Inserting 1 = 1+ + 1~ between p m i and B 
in (kem\p m > B\lem) and between B and p m > in 
(lem\Bp m i \kem) gives 

(kem\p m iB\lem) = (kem\lmPm,'B\lem) 
+ (fcem|l~/w B\lem) 



(lem\Bp m '\kem) = ^em\B\^ n p m i\kem) 
+ (Zem|i?l~/5 m / \kem) 



and 



In Sec. IC 21 we show that the second terms in the above 
expressions vanish, i.e. 

(kem\l-p m ,B\lem) = (C12) 
(lem\Bl-p m ,\kem) = 0. (C13) 

On using 1+ = J2k'e' \k'e'm)(k'e'm\ from Eq.|C4| the 
terms involving 1+ are evaluated as 



Note that these expressions arc finite only for m' > m. 
Using the definition of the reduced density matrb*^ 



m -4m / 
Pied 



l (k,k') = Tr e [(kem\p m '\k'em)], 
we arrive at the following expression for G™/ (lu) 

Pred l K ' K ) D k'l 



m —1 

E J2 A 

m-mo+l Ikk 



Ik' 



LU ■ 



Ef 



EF + id 



m — 1 

E E^s 

m=m a + l Ikk' 



Pred \ K ; K ) a lk' 

u + E™ -E™ + iS 



Hence, the final expression for G(lu) = 
E5= mo+ i ™rn>(GU") + Gt,{u) + G*f,(w)) is given by 



G(lu) 



m IV 



+ 



+ 



+ 



+ 



N 

E 

m'— mo + 1 
7V-1 

E f^E^i 

m'-mo + l Ik 
N-l 

W, 



'hi 



E 



T^Y A klB?k 
m'-mo+l m kl 
N-l 

E E^; 

m=m + l Ikk' 



(e- 




+ 


' ) 


LU-\ 


-Ef 


- 4 


-id 




e~ 






LU-\ 


-Ef 


- Ef + 


■id 




e~ 







El 



Ef 



id 



N-l 

E 



m=m + l kk'l 



L—f U LU- 



Rred(^fc , )^["/ 

+ Ef 1 -E™ + id 

R-redV^ J k)E"k' 



Ef 



E™ + id 



where, in the last two terms we rearranged the summa- 
tions over m' and m and introduced the full reduced den- 
sity matrix 



R-rod(^) & ) 



N 

E 

rn'— m+ 1 



m — >m 
w m' Fred 



(k,k'). 



(kem\lmf>m> B\lem) = ~y^^(kem\p m i\k' e'm) (k' e'm\B\lem) 

k' e' 

= y y^(kem\p m '\k' em) (k' em\B\lem) 

k' 

= ^2(kem\p m ,\k'em)B]!} l 

w 

(lemlBl^pmi \keni) = ~^^(lem\B\k / e'm) (k'e'm\p m t \kem) 

k' e' 

= (lem | B | k 1 em) (k'em \ p m i \ kern) 

k' 

= yB%,(k'em\p m ,\kem) 



Note that the meaning of this quantity is completely anal- 
ogous to the reduced density matrix introduced by Hof- 
stetter in Ref. [36| except that one obtains reduced den- 
sity matrices at iteration m by eliminating environment 
degrees of freedom e = e m j r \e m j r i---&N from the FDM 
(|C7[) instead of the density matrix for iteration N. In 
addition, the former is built from the complete set of 
eliminated states, as opposed to the retained states of 
iteration N in the approach of Ref. HH. The above ex- 
pression for G(lu) is identical to that in Ref. [H. We 
have checked that the sum rule for the spectral function 
A„{oj,T) = -±Im[G AB (u))} 



A a (ui, T)duj = 1 



is satisfied exactly (to machine precision) when using the 
discrete (unbroadened) form of the spectral function as 
in Ref. ^3. 



19 



1. Calculation of weights «;„ 



The expression for w rn in (|C11|) involves Z which con- 
tains eigenvalues from all iterations ml = m + 1, . . . , N. 
In evaluating these expressions, one should therefore use 
the absolute energies for the E™. Since, in practice, the 
iterative diagonalization of the Hamiltonian H m involves 
subtraction of g roundstate energies and rescaling at each 
m (see Ref. l30f ) , one has to keep track of the subtracted 
groundstate energies and return to the actual physical 
energies relative to a common absolute energy reference 
in evaluating w m and Z. We take this absoute energy ref- 
erence to be the ground state energy of the last Wilson 
iteration N. Thus, if E^ s is the true groundstate energy 



of H m , we use E\ n -> E™ 
evaluating w m /Z m and Z 



GS 



Z n 



-? E csZ m in 



Z n 



z 



T N , . , 4 N - m 'e-P E GsZ m , 

L^m' — mo + 1 m 
m'— mo + 1 



2. Proof of Eq. (fC12j) and Eq. (fC13l) 

Using the expression for p m i we easily find that 



E E 

ra" — mo + 1 le 



lem"){lem"\ x 



Y\l' e ' m ')^ (/Vm'l 



l'e 



E 



J m m 

\ p m > if ml < m; 
10 if m! > m . 



E \ lem '} 



Z,i 



-(lem'\ 
(C14) 



Hence (fcem|l~p m /i?|^em) in Eq. (|C12j) involves matrix 
elements of the form (kem\l'e'm') for m! < m, which van- 
ish, since all retained state at iteration m have no overlap 
with eliminated states at iterations m! < m (i.e., elimi- 
nated states of previous iterations are not used to obtain 
retained states of later iterations) . The same arguments 
can be used to prove Eq. (|C13[) . 



Appendix D: Thermal conductance and 
thermopower of quantum dots 

For completeness, we outline here the derivation of 
thermoelectric transport through a strongly interacting 
quantum do t 21 ' 22 . The electrical, II, and heat cur- 
rent, I®, from the left lead to the quantum dot can 
be expressed in terms of the particle number Nl = 



Efc<x C Lka C Lka and energy H L = £ fe(J ^Lkac{ ka c Lk(T of 

the left lead, via 



II 
I" 



-eN T 



(Dl) 



L = H L - ii L N L = — \H L - vlN l ,H], (D2) 

where H is the Hamiltonian ([TJ. In terms of the lesser 
Green's function's G^ kLa (t,t') = i{c{ La (t')d a (t)) and 

= i(dUt')c k L*(t)) = -(Gl, kLa y, the 
above currents are given by 



2e 



Re 



E^ 



fcLi 



■(*,*) 



da.kL 



(D3) 



.(D4) 



The lesser Green's function G da kLa (t,t') can be ex- 
pressed via equations of motion solely in terms of Green's 
functions of the dot and the non-interacting Green's func- 
tion for the left lead. After some lengthy algebr a 22 ' 40 , 
one finds the following expressions for the currents in 
terms of the retarded, G r dt7 = Gda(uJ + iS), advanced, 
G da = Gd a {u) — iS) and lesser Green's function, G da {uj) 
of the dot: 



II = 



dcof L [(G<» 

G d a - G a da )] , (D5) 
dw(u> - (j, L )f L [G^ a (uj) 



f L (u)(G d(7 -G a d<7 )] 



(D6) 



where f L (w) = (l+e-^"^)/^^))" 1 is the Fermi func- 
tion of the left lead and Tl = nNptW is the hybridiza- 
tion strength of the dot to the left lead as defined in 
Sec. UU By using current conservation II = —Ir, one 
can eliminate the lesser Green's function from the above 
expressions to arrive at the final expressions used in this 
paper 

II = |E/ Mh-fR)TM, (D7) 
1 l = -JE/ ^(^-^)(/i-/fl)T rfCT H.(D8) 

The quantity Tda(u) acts as a transmission function and 
is given by 

T d » = 2iJ^-(G da -G a da ) (D9) 

1 L + 1 R 

The electric and heat currents are expanded to linear 
order in ST = T L - T R and SV = V L - V R 



h ) f in Lis \ [ 6V 



(D10) 
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defining, thereby, the transport coefficients Lij,i,j = 
1,2. In terms of the latter, the transport properties are 
given by 



G(T) 
S(T) 

K C {T) 



lim Il/5V\ S t=o = L n 

SV^>0 



= - lim SV/8T\ Il=0 

<5T->0 

= - lim I Q /6T\ lL=a 

8T->0 

= L12L21/ Ln — L 22 . 



(Dll) 
L 12 /L n (D12) 

(D13) 



Finally, the Lij are simply expressed in terms of the fol- 
lowing transport integrals 



In(T) = lj 'de e n T(e) 



df 
de 



(D14) 



via in = g 2 Iq,L2i = —elx/T and L22 = I2/T. Substi- 
tuting these values for L.y into (|D1HID13|) results in the 
expressions given in the text. 
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